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Introduction 

The question of the existence of thermodynamic hmit for all standard (i.e, 
short range, see Remark 2 below) models of spin glasses with two-body in- 
teractions has been settled long ago by Khanin and Sinai |KSj and later 
generalized to more general interactions by Zegarlinski [Zej (previous refer- 
ences on the subject include jLe], |Vuj - jPFj ). The sharper property of the 
monotonicity of the free energy in the volume has been proved by van Enter 
and van Hemmen |EHj . In the long-range case (the most important examples 
being the Sherrington-Kirckpatrick model |MPVj , the REM and the GREM 
|DGp the existence and the monotonicity has been instead proved only very 
recently |G1] . [nT2] . |CDGGj . The proof relies on an interpolation argument 
introduced in |GTj which has the advantage of yielding the subadditivity of 
the free energy (equivalently, superadditivity of the pressure). Exactly as 
in the ferromagnetic case |Ruj . and in |EHj . the subadditivity entails the 
important property of the monotonicity of the free energy (pressure) as the 
volume increases. 

In this paper we show that the above interpolation argument can be 
applied (actually in a slightly different form) to the short range case. For 
Gaussian couplings, in the summable case and in the non-summable one as 
well, we generalize the Khanin-Sinai and van Enter-van Hemmen results. 
Namely, for general, mixed, short-range n— body interactions (n arbitrary: 
for the physical relevance of > 2 see e.g.jLi]), and assuming free bound- 
ary conditions, the free energy, internal energy, ground state energy are not 
only bounded but also decreasing in the volume. Hence their thermodynamic 
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limit is reached monotonically. We remark that in the disordered case the 
monotonicity is even more relevant than in the ferromagnetic one because 
the ground state energy is tacitly assumed monotonic in all numerical simu- 
lations; for a discussion of this point, see e.g. [Ri], |BCDGj . 

The conclusion which may be drawn by this paper, together with referen- 
ces jGTj and |CDGGj . is that as far as the thermodynamic limit is concerned 
Gaussian spin glasses in free boundary conditions do not differ from ordinary 
ferromagnets: in both cases pressure, internal energy and ground state en- 
ergy are bounded and monotonic in the volume. 

Definitions and Examples 

Let M be a countable set and consider a finite subset A C M of cardi- 
nality |A| = A^. To each element i G A we associate a dynamical variable 
CTj G iS C (for some fixed integer k) equipped with an a priori proba- 
bility measure z/j. For each X C A we consider ax = {o"i}jex and a function 

In analogy to |Ruj (Sect. 2. 4, formula 4.3) and [Zej we define the random 
potential as 

f/A(J,a)= Y^Jx'^xM , (1) 

XcA 

(with $0 = 0) under the following assumption: the coefficients Jx are inde- 
pendent Gaussian variables with zero mean and variance depending only on 
X (and not on A) 

Av( Jx) = , Av(4) = . (2) 
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Examples: 



Here M = Z'', and A is a cube. 

1. The Edwards-Anderson model. S = {+1, -1}, iy{ai) = The 

nearest neighbor case is defined by $Ti,n'(^n) <^n') = <^n<^'n ^oi \n—n'\ = 1, 
$x = otherwise, and A'j^ = (? . More generally one may consider a 
short range interaction with with — \n — n'l"^''", a > 1/2, or a 
many-body interaction with a suitable decay law. 

2. Multicomponent spin models (Potts models): S — {1, 2, q}, i^(crj) — 



equal and zero otherwise. 
3. Continuous spin models: S — M/^, v{ai) — dix{x) > 0, / dii{x) ~ 1 



(unbounded case) or 5 = T'^, t^((Tj) = dcp (bounded case); 
4. Lattice gases: here S — {0, 1}, v{ai) = ^[Sq + Si]. 
Remcirks: 

1. Of course the examples may be considered on every finite dimensional 
lattice hke "L^ or the triangular lattice etc. 

2. The property that is independent of the volume A characterizes 
the short range case, such as the Edwards- Anderson one. In mean 
field (long range) models, such as the Sherrington-Kirckpatrick one, 
the variance has to decrease with N in order to have finite energy 
density. 




1=1 



8i , $x(cx) = S^x where 5^^^ = 1 if all components of ax are 
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Denoting Pj\^{da) — Yli^A^i^ii'^i) we define: 

1. The random partition function 

Za{J) := / PA(da)e^-(^''^) , (3) 

2. The random Gibbs-Boltzmann state 

f Pa (da) - e'^AC^.o-) 

3. Tfie quenched state 

< - > Av(a;(-)) , (5) 

4. The quenched pressure 

Pa:= Av(lnZA(J)) . (6) 

5. The quenched potential 

Uj,:^<UA{J,a)> (7) 

We remind that the free energy Pa is — /3~^Pa, and the internal energy Pa 
is 

SUPERADDITIVITY 

Lemma 1 

< Jx^x > > . (8) 
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Proof. 

We remind the integration by parts for Gaussian variables 

Av(J^/(J)) = A|Av(|i^) , (9) 

and the correlation derivative formula 

'^^^=ujm-u.{^xf >0. (10) 

By applying successively © and (fTUj) we obtain 

< Jx$x > = Av(Jxcj($x)) = 

= A^Av(^($^)-^($x)') >0. (11) 

As a corollary of lemma 1 we have 

< [/a( J, a) > = ^ A^Av {uj{^\) - uj{^xf) > (12) 

XcA 

Definition 1 Consider a partition of A into n non empty disjoint sets Ag/ 

n 

A=UA,, (13) 

s=l 

A,nA,, =0. (14) 

For each partition the potential generated by all interactions among different 
subsets is defined as 

n 

Ua = Ua-J2u^^; (15) 

s=l 

from (OJj we have that 

Ua=Y1 -J^^x (16) 
where Ca is the set of all X G A which are not subsets of any Ag. 
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Theorem 1 The quenched potential is supemdditive: 

n 

< f/A > > 5^ < f/A. > (17) 
s=l 

Proof. Direct consequence of (jH)). In fact: 

< f/^ > = ^ < J^$^ >= ^^^^ (^(*^) - ^^^^f) ^ ■ (18) 
xeCA X£Ca 

Theorem 2 The quenched pressure is superadditive: 

n 

Pa>$^Pa. . (19) 

Proof. 

To each partition of A we associate the interpolating potential for < t < 1 

n 

u^it) = Y,vt:ui'], (20) 

s=0 

with to = t, ts = {1 — t) for 1 < s < n, U^^ = U\ and 

< = E Jx^^ ' (21) 

XcAs 

(s) 

where any J)^ is a centered independent Gaussian 

Ay(^jP4'^)=6sJx,y^x (22) 

(the symbol Av is here the average with respect to all the J's). We define 
the interpolating partition function 

ZA(t) = / PA(rfa)e^-W , (23) 
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and we observe that 

n 

Za(0) = n^A.(/^)) , ^a(I) = Z^{J) . (24) 

s=l 

Consider the interpolating pressure 

PaW := Av(lnZA(t)) , (25) 
and the corresponding states 0t{—) and < — >t. Thanks to we get 

n 

Pa(0) = 5^Pa. , Pa(1) = Pa. (26) 



s=l 

We observe now that 



with eo = 1 and = — 1 for 1 < s < n. For each s we have 

<uS>t= E <Jx^x>t; (28) 

XcAs 

Integrating by parts each addend we obtain 

< 4^)$x >t = v^A^Av (^i(<l>^) - ^t($x)') (29) 
and summing up all the contributions in (|27|): 

|pA(t) = J]A|Av(c.,($^)-a;,(<|.x)^)- 



5^ 5^ A^Av(.;,(<|.^)-.;,($;,)2) 

s=l XcAs 

J2 AxAv (a;t(<l>|) - ujti^xf) > . (30) 



xeCA 



IFtoiu. (j26p and (|3U|) we immediately get formula (jl9|) . 



BOUNDEDNESS 

For any random potential we define the quantity 



I \U\ I = sup i-Av {UaU a)') = sup ^ 5^ A^$^ . (31) 



XcA 



Potentials with a finite ||?7|| are called stable. 

Theorem 3 A stable random potential admits an internal energy and a 
quenched pressure bounded by the volume. 

Proof. By 

<UA{J,a)>= J2^x^''H'^x)-^{^x)^) <2\\U\\N. (32) 

XcA 

Using the Jensen inequality 

Pa = Av(lnZA(J)) < lnAv(ZA(J)) = 
= \n [ PA(rffT)Av (e^'^W")) = \n [ P^{da)e-^^^cA^'x^'x = 

< l\\U\\N (33) 
As a consequence for finite \ \U\ \ one has 

sup ^Ua < oo , (34) 

A iV 

and 

sup ^Pa < oo . (35) 

A iV 

Thermodynamic Limit 

Let us verify the stability condition in the above examples. 
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1. Edwards- Anderson. 



For the nearest neighbor case 



J2 ^x = 2dNC^ 



(36) 



(n,n') 



2. More generally for the short range case with a > 1/2 



1 



(37) 



l\2da — 



< const N 



n — n' 



By Theorems 1 and 2 the previous models have an internal energy per par- 
ticle and a free energy per particle which exist in the thermodynamic limit. 

Remarks: 

1. We point out that for short range models we only need to impose 
the boundedness condition (j31|) while the superadditivity always holds 
thanks to the condition of independence of the variance from the 
volume (see also the remark 2). In the mean field case the variance of 
the interactions depends on the volume and subadditivity is based on 
an inequality among the covariances 



One may check that such an inequality reduces, in the short range case, 
to the positivity of the right hand side of (jnOJ. 

2. Our result may be extended in two directions by exactly the same pro- 
cedure of jGTj . First one can prove by standard probability arguments 



NiCN,{a, r) + N2CN,{(r, r) - NcN{(r, r) > . 



(38) 
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that the above statement entails the almost sure convergence of pres- 
sure and ground state energy per particle. Second our result may be 
extended to non Gaussian J (see section 4.2 of |GTj and jTj): if Jx 
is for all X an even random variable with a finite 4th moment the 
integration by parts Q is replaced by the more general formula 

Av(JxF(J)) = Av(4F'(J)) - \m\Jx\ (4 - x^)F"\x)dx) . 

4 J-\Jx\ 

(39) 

When used within Theorem 2 it generates in formula (j3(Jj) a correction 
of order 0{^/N). Once the density is taken the correction vanishes in 
the thermodynamic limit. In general, however, we cannot establish its 
sign and the monotonicity is lost. 

3. The present result holds for free boundary conditions. In general 
it is proved [Ze] by standard surface over volume arguments that the 
quenched quantities are independent of the boundary conditions, but 
the monotonicity property is lost. 

4. It is interesting to observe that the interpolating strategy does ap- 
ply also to standard ferromagnetic systems. Consider for instance the 
d-dimensional ferromagnetic Ising model with nearest neighbor Hamil- 
tonian H^{a) = —J2{nn')^riC''n'- An interpolating functional would 
be 

a{t) = log^e-'^[*^^M+(i-*)S"-^^=('^)l . (40) 
An easy calculation which goes parallel to Theorem 2 yields 

^(t) = Yl ■ (41) 

(n,n')eC 
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Since < t < 1 the t-interaction on (j4(Jj) is still ferromagnetic and 
the Griffiths inequality (see for instance [Ruj) 0t(o"nCr„/) > gives the 
positive sign of the former expression ensuring the monotonicity of the 
limit. 
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